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FOREWORD 


This  document  was  prepared  for  the  Special  Communications  Project 
Office  of  the  U.S.  Naval  Electronic  Systems  Command  by  IIT  Research  Insti 
tute  under  Contract  N00039-76-C-0141. 

The  technical  effort  reported  herein  is  intended  to  support  bio- 
logical studies  in  connection  with  the  Navy's  Seafarer  ELF  communications 
system.  The  report  covers  an  analysis  of  the  electric  field  levels  in- 
duced in  individual  tissue  culture  samples  when  exposed  to  uniform  elec- 
tric and  magnetic  fields  in  a laboratory  simulator.  Calculations  are 
included  to  derive  sample  field  levels  for  the  NMRI  tissue  culture  experi 
ment.  A companion  report  "ELF  Electric  and  Magnetic  Field  Simulation  for 
a Laboratory  Biological  Experiment,"  IITRI  Technical  Report  E6357-10,  de- 
scribes the  field  simulator  developed  by  IITRI  for  the  NMRI  experiment. 
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1.  INTRODUCTION 

Modern  technology  has  developed  many  applications  for  radio  frequency 
and  microwave  electromagnetic  energy.  These  applications  include  communi- 
cations, remote  detection,  industrial  processing,  food  preparation,  and 
medical  diagnosis  and  therapy.  This  extensive  usage  has  led  to  concern 
about  the  potential  for  biological  hazards  from  electromagnetic  radiation 
exposure  and  to  the  setting  of  national  standards  for  the  maximum  permis- 
sible exposure  for  humans. * The  possible  effects  on  organisms  of  extremely 
low  frequency  (ELF,  20  to  100  Hz)  electric  and  magnetic  fields  are  also 
being  studied.  This  research  includes  the  Navy's  proposed  ELF  Communica- 
tions System^  as  well  as  power  systems. ^ The  research  spans  a broad 
spectrum  of  biological  sciences,  from  enzymes  to  mammalian  behavior,  and 
from  bacteria  and  slime  molds  to  primates. 

Recently,  research  has  been  initiated  at  the  Naval  Medical  Research 
Institute  (NMRI ) to  study  the  possible  effects  of  the  electric  and  magnetic 
fields  associated  with  the  Navy's  proposed  ELF  Communications  System 
(Seafarer)  on  human  tissue  cultures.  The  necessary  electrical  engineering 
support  for  this  experiment  is  provided  by  the  I IT  Research  Institute 
(IITRI).  The  apparatus  used  to  generate  and  control  the  required  electric 
and  magnetic  fields  was  designed,  fabricated  and  installed  by  IITRI.  In 
support  of  this  design,  an  analysis  was  performed  to  determine  the  field 
internal  to  the  tissue  culture  as  a function  of  applied  field  and  exposure 
geometry.  The  purpose  of  this  report  is  to  present  the  details  of  this 

5 

analysis.  A companion  report  discusses  the  design  and  operation  of  the 
electric  and  magnetic  field  simulator. 

The  ELF  electromagnetic  field  simulator  designed  for  tissue  experi- 
ments at  NMRI  employs  parallel  plates  to  generate  a uniform  electric  field 
and  a multi -coil  to  provide  a uniform  magnetic  field.  A photograph  of 
this  simulator  is  shown  in  Figure  1.  The  tissue  under  test  is  normally 
placed  in  a plastic  tissue  culture  tray.  The  photograph  of  Figure  2 
illustrates  one  of  these  tissue  culture  trays. 

In  the  test  volume  of  the  simulator,  both  electric  and  magnetic 
fields  are  induced  inside  the  biological  tissue.  The  important  aspect 
of  the  present  problem  of  electromagnetic  coupling  and  interaction  is 
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Fig.  2 A TISSUE  CULTURE  TRAY  FOR  THE 
EXPERIMENTS  AT  NMRI 
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the  determination  of  the  electric  field  strength  and  the  magnetic  field 
strength  induced  inside  the  tissue  culture.  For  the  case  of  present 
interest,  the  magnetic  field  inside  the  tissue  culture  has  a magnitude 
equal  to  that  of  the  impressed  field,  and  suffers  no  significant  distor- 
tion because  the  frequency  of  interest  is  low  and  the  permeabilities  of 
the  biological  tissue  and  their  surrounding  medium  are  essentially  that 
of  free  space.  To  determine  the  electric  field  strength  internal  to  the 
biological  tissue,  we  must  either  conduct  an  experimental  measurement  or 
develop  a theoretical  analysis.  For  the  periodic  lattice  structure  (see 
Figure  2)  considered  here,  it  is  difficult  to  accurately  measure  the 
internal  electric  field  strength  because  of  the  fact  th3t  the  electric 
field  strength  internal  to  a very  high  permittivity  object  is  low,  and 
because  the  test  specimens  are  small.  As  a result,  the  electric  field 
coupled  into  the  biological  tissue  in  the  plastic  tray  of  Figure  2 must 
be  determined  analytically. 

The  purpose  of  this  study  then  is  to  analytically  determine  the 
electric  field  induced  inside  biological  objects  which  are  embedded  in- 
side a dielectric  slab  in  a periodic  fashion.  Specifically,  we  want  to 
know  quantitatively  the  electric  field  strength  inside  the  tissue  culture 
when  the  tissue  culture  tray  is  exposed  to  the  uniform  electric  field 
generated  by  the  simulator. 

The  biological  object,  characterized  by  the  constitutive  parameters 
ej,  Up  is  assumed  to  be  homogeneous  and  isotropic  where  is  the  permit- 
tivity and  Uj  is  the  permeability  of  the  biological  tissue.  The  surround- 
ing dielectric  medium  in  which  the  objects  are  embedded  is  also  homogeneous 
and  isotropic  with  parameters  Ug-  either  medium,  the  permittivities 
could  be  arbitrary  and  complex,  and  the  permeabilities  are  considered  to  be 
that  of  free  space:  u^  = U£  = PQ* 

Two  analytical  methods  are  developed  for  the  solutions  of  the  electric 
field  induced  internal  to  the  tissue.  The  first  method  is  based  on  the 
Lorentz  theory^  developed  in  the  analysis  of  an  artificial  dielectric.  Anal- 
ytical expressions  are  obtained  for  the  induced  internal  electric  field.  The 
second  method  uses  an  equivalent  circuit  approach.  While  the  first  method 
is  more  accurate,  the  solution  of  the  second  method,  although  admittedly 
approximate,  is  simpler  to  obtain  and  easier  to  understand.  The  numerical 
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calculations  show  good  agreement  between  the  results  of  these  two  methods. 
The  calculated  results  for  the  internal  electric  field  are  graphically 
presented  for  certain  lattice  structures  of  the  culture  tray  and  for 
various  geometric  shapes  of  the  object. 

For  the  most  part,  the  details  of  theoretical  analysis  are  presented 
in  the  appendices,  whereas  the  important  results  are  summarized  in  the  text 
of  this  report. 
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ANALYTICAL  METHODS  AND  SOLUTIONS 


As  shown  in  Figure  1,  the  simulator  for  use  in  tissue  culture  exper- 
imentation employs  parallel  plates  to  generate  an  electric  field,  and  a 
multi-coil  to  generate  a magnetic  field.  The  basic  design  requirement  of 
this  simulator  is  that  it  provides  a uniform  electromagnetic  field  over  a 
desired  test  volume.  A design  goal  was  that  variation  of  the  field  within 
the  volume  should  not  be  more  than  10  percent.  Therefore,  the  impressed 
field  considered  in  this  study  is  assumed  to  be  uniform.  In  addition, 
only  the  electric  field  coupling  is  considered,  because  the  permeabilities 
of  the  biological  objects  and  their  surrounding  medium  are  considered  to 
be  the  same.  Consequently,  the  magnetic  field  induced  inside  the  object 
suffers  no  distortion  or  variation  from  the  impressed  magnetic  field. 

The  electric  field  inside  the  object  is  induced  by  the  impressed 
electric  field.  The  tissue  cultures  which  are  embedded  in  the  tissue 
culture  tray  for  testing  are  exposed  to  the  uniform  electric  field  gener- 
ated by  the  simulator.  A sketch  of  the  geometry  is  given  in  Figure  3. 
oecause  of  the  periodic  arrangement  of  the  tissue  culture  tray,  the  cou- 
pling mechanism  is  more  complicated  than  that  of  a single  isolated  object. 
To  determine  the  electric  field  strength  internal  to  the  object,  we  must 
take  into  account  the  interaction  of  the  fields  among  all  the  objects. 

The  solutions  can  be  obtained  from  the  interaction  field  approach  pre- 
sented below.  A simple  approximate  method  based  on  equivalent  circuit 
technique  is  also  given. 

2.1  Interaction  Field  Approach-Lorentz  Theory 

When  the  interaction  of  a uniform  electric  field  with  a single  iso- 
lated object  is  considered  at  extremely  low  frequencies,  the  coupling 
mechanism  is  essentially  governed  by  the  field  boundary  conditions.  The 
coupling  is  stronger  for  objects  having  more  elongated  geometries.  The 
biological  tissue  culture  tray  of  Figure  3 considered  here,  however,  is  a 
large  scale  model  of  an  actual  dielectric  obtained  by  arranging  a number 
of  identical  biological  tissues  (objects)  in  a regular  three-dimensional 
(or  two-dimensional)  pattern.  Thus,  the  coupling  mechanism  becomes  more 
complicated  resulting  from  the  interaction  of  the  fields  among  all  the 
objects.  To  take  the  interaction  field  into  account,  we  note  the  Lorentz 


theory  which  was  originally  developed  for  the  theory  of  a dielectric. 

This  theory,  which  has  been  extensively  used  in  the  study  of  artificial 
dielectrics,  considers  dipole  interactions  among  objects,  and  provides 
accurate  results  for  object  spacings  limited  to  less  than  0.1  of  a wave- 
length and  for  object  sizes  small  compared  to  a wavelength. 

The  electric  field  coupling  to  biological  objects  can  be  described, 

O 

as  Stratton  did,  in  terms  of  polarization  (or  dipole  moment)  because  when 
exposed  to  a uniform  electric  field,  electric  charges  are  induced  inside 
the  object.  These  electric  charges  are  displaced  inside  the  objects  so  as 
to  set  up  induced  fields  which  will  eventually  satisfy  the  boundary  condi- 
tions at  the  object  surface.  Since  an  electric  dipole  or  simply  a dipole 
is  a name  given  to  two  point  charges  of  equal  magnitude  but  opposite  sign 
separated  by  a small  distance,  the  term  polarization  (or  dipole  moment) 
can  be  used  to  characterize  the  nature  of  biological  objects  under  the 
influence  of  an  impressed  electric  field.  The  ELF  electric  field  coupling 
to  dielectric  objects  of  various  geometric  shapes  is  discussed  in  terms  of 
polarization  in  Appendix  A. 

For  the  present  problem,  it  is  convenient  to  consider,  without  loss 
of  any  generality,  that  the  biological  tissues  have  their  principal  axes 
of  polarization  coincident  with  the  coordinate  axes  (as  shown  in  Figure  T) 
so  that  the  induced  dipole  moment,  p,  of  the  object  (tissue)  is  in  the 
same  direction  as  that  of  the  impressed  field.  At  extremely  low  frequencies, 
the  field  acting  to  polarize  any  object  may  be  assumed  to  be  a uniform  field 
equal  to  the  field  at  the  center  of  the  object.  When  each  (biological)  ob- 
ject is  replaced  by  its  equivalent  dipole  moment,  p,  (obtained  from  Appendix 
A),  the  tissue  culture  tray  of  Figure  3 considered  here  reduces  to  an  equiv- 
alent dipole  array.  This  equivalent  dipole  array  is  shown  in  Figure  4. 

Since  each  object  exhibits  a dipole  field,  interaction  occurs  between 

objects.  This  interaction  effect  may  be  evaluated  in  terms  of  a constant, 

9 

B,  called  the  interaction  constant  which  has  been  defined  by  Collin  ’ in 
the  work  in  connection  with  artificial  dielectrics  and  aperture  coupling 
in  waveguides.  The  interaction  constant  for  the  present  problem  has  been 
derived  in  Appendix  B using  the  image  theory  for  potential  functions  of 
either  a three-dimensional  array  or  a two-dimensional  array. 
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Fig.  4 AN  EQUIVALENT  DIPOLE  ARRAY  FOR  THE  PERIODIC 
LATTICE  OF  FIGURE  3. 
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As  a result  of  this  interaction,  the  net  combined  effect  of  all  the 
objects  is  to  produce  a new  average  dipole  moment  at  the  center  of  each 
object.  With  the  polarizabil ity  constant  (Appendix  A)  and  the  interaction 
constant  (Appendix  B)  obtained,  the  induced  electric  field,  , internal 
to  the  object  under  the  influence  of  a uniform  electric  field  may  be  ob- 
tained by  using  the  concept  of  polarization  and  the  theory  of  artificial 
dielectrics.  The  results  obtained  in  Appendix  C are 


$ae/(C-l)  } 

Ei  = ^(l-aeB)kv^  Ep 


(1) 


where 


U2+  (£e/£o)Ul  + V} 


(2) 


e 


e 


ag/abc  J 


(3) 


In  these  expressions,  the  following  parameters  and  notations  are  used. 

a,b,c:  object  spacing  along  the  direction  of  x,y,z  axis, 
respectively,  m. 

£,(£3):  the  spacing  between  the  tissue  culture  tray  and  the 
upper  (lower)  plate  of  the  simulator,  m. 

£ 2 : thickness  of  the  periodic  lattice,  m. 

V : voltage  applied  between  the  parallel  plates  of  the 
simulator,  volts. 

-12 

eQ:  free  space  permittivity,  8.854  x 10  farads/m. 

e,:  permittivity  of  the  biological  tissue,  farads/m 

£p:  permittivity  of  the  medium  outside  the  biological 
tissue,  farads/m. 


S 


Ve2- 

volume  of  the  biological  tissue,  m . 


ot 


e 


polarizabil ity  constant  (the  polarizability  constants  for 
objects  of  various  geometric  shapes  are  listed  in  Table  A-2 
in  Appendix  A). 
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g:  interaction  constant  (the  interaction  constants  for  the 
three-dimensional  dipole  array  and  the  two-dimensional 
array  are  given  in  Equations  (B-7)  and  (B-ll)  of  Appendix  B, 
respectively) . 

Using  »e  and  & obtained  in  Appendices  A and  B,  Equations  (l)-(3)  can 
be  used  to  calculate  the  induced  electric  field  internal  to  the  biological 
tissue,  E^ . As  an  illustrative  example,  let  us  consider  the  biological 
objects  having  a prolate  spheroid  geometric  shape.  Using  the  expression 
of  ag  listed  in  Table  A-2  of  Appendix  A,  it  can  be  easily  shown  that 
Equations  (1)  and  (3)  become 


i " {l-kveU-l)/A} 

h + kv(5-l)(abc-6)/A| 
ee  = e2j  l-kveU-l)/A  \ 

respectively,  where 

A = {l+fe-l)(u2-l)[ucoth'1  (u)-l] } 


(4) 


(5) 


(6) 


u 


-1/2 


(7) 


★ 

In  Equations  (4)-(7),  r and  r.  denote  the  major  and  minor  radius  of  the 

a D 

prolate  spheroid,  respectively,  and  Ep  again  is  given  by  Equation  (2). 

This  result  indicates  than  when  Ej  = C=l)>  we  have  eg  = 

and  A = 1 which,  in  turn  gives  E,.  = V /[j^  + (eg/Cg)  • This  result 

is  identical  to  that  obtained  for  a homogeneous  dielectric  slab  situated 
between  a parallel  plate  with  an  applied  voltage  V . It  is  also  interest- 
ing to  note  that  when  rb  decreases,  A decreases  accordingly  so  that  when 

rb  = 0,  the  result  of  Ei  reduces  to  the  result  for  a very  thin  and  long 

cylindrical  object.  For  r./r.  = 1,  we  find  A approaches  1/3  and,  hence, 

the  result  of  E^  becomes  that  of  a biological  sphere. 


★ 

In  Appendix  A,  the  major  and  minor  axes  of  the  prolate  (or  oblate)  spheroid 
are  designated  by  a and  b,  respectively. 
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2.2  Circuit  Approach 

In  this  section,  a simple  circuit  approach  is  discussed.  While  this 
method  provides  only  approximate  results,  the  solutions  are  useful  for  a 
quick  estimate  of  the  electric  field  induced  inside  the  biological  tissue 
and  provide  for  an  understanding  of  the  coupling  phenomenon. 


As  noted  previously,  the  permittivities  of  the  biological  tissue 
and  its  surrounding  medium  could  be  complex.  Denoting  = e.'-jej" 
and  ^ = e2  ^e2~*  t^ie  condl,ct1 vities  of  the  tissue  and  the  medium  are 
a ^ and  02  = ■>  respectively.  Hence,  the  real  part  of  the 

permittivity  (e^  or  Eg"),  like  a capacitor  represents  the  stored  energy, 
and  the  imaginary  part  of  the  permittivity  or  e 2")  is  responsible 

for  the  power  loss.  It  can  easily  be  seen  that  the  tissue  culture  tray 
and  the  parallel  plate  simulator,  shown  in  Figure  3,  may  be  represented 
by  an  equivalent  circuit  obtained  from  an  approximate  geometric  model. 

Two  realistic  geometric  models  which  approximate  the  tissue  culture  tray 
of  Figure  3 are  illustrated  in  Figure  5,  along  with  their  corresponding 
equivalent  circuits.  The  two  geometric  models  presented  here  are  consid- 
ered to  be  illustrative  and  provide  reasonably  good  approximations. 


In  the  tissue  culture  tray,  the  medium  which  surrounds  the  biological 
tissue  normally  is  a thin  layer  of  perfect  dielectric  (i.e.,  e2~  = 0). 
Hence,  the  capacitor,  Co,  is  sufficient  to  characterize  this  medium.  The 
same  is  true  for  the  free  space  capacitor,  C^.  The  capacitance  and  resis- 
tance in  the  equivalent  circuits  can  be  obtained  from 
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In  Equation  (8),  we  note  that  SR  + S^.  = 1,  where  SR  is  the  ratio  of 
the  cross-sectional  area  of  the  tissue  to  that  of  the  tissue  culture  tray 
in  the  x-z  plane  of  Figure  5,  and  Sc  denotes  the  ratio  of  the  cross-section- 
al area  occupied  by  the  plastic  (dielectric)  to  that  of  the  tissue  culture 
tray. 

From  these  equivalent  circuits,  the  voltage  across  the  resistor  of 
either  (in  Figure  5a)  or  R (in  Figure  5b)  can  be  obtained  easily.  If 
this  voltage  is  VR,  then  the  electric  field  induced  inside  the  biological 
tissue  is  E..  = On  this  basis,  the  equivalent  circuit  of  Figure  5b 

is  expected  to  yield  a lower  bound  estimate  of  the  field,  Ei , since  the 
flux  would  be  uniform  in  this  model.  The  equivalent  circuit  of  Figure  5a 
will  yield  an  upper  bound  estimate  of  E^  because  it  assumes  that  the  flux 
captured  by  the  tissue  depends  only  on  the  cross-sectional  area  ratio. 
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3. 


CALCULATED  RESULTS 


Although  the  analysis  presented  in  this  study  is  generally  applicable 
for  calculating  the  induced  electric  field  inside  the  biological  tissue  for 
any  dimension  of  the  simulator  and  for  objects  of  any  geometric  shape,  the 
calculations  presented  here  emphasize  the  results  obtained  for  the  specific 
simulator  and  the  tissue  culture  tray  fabricated  for  the  NMRI  experiment. 
More  parametric  calculations  are  presented  in  Appendix  C. 

The  simulator  considered  here  generates  ELF  fields  at  a center  fre- 
quency of  76  Hz.  The  complex  permittivity  of  the  biological  tissue  is 
considered  to  that  of  a high  water  content  tissue.  A search  of  available 
literature  did  not  reveal  an  accurate  value  of  the  permittivity  for  the 
high  water  content  biological  tissue  at  extremely  low  frequencies.  For 
higher  frequencies  above  1 kHz,  the  value  of  permittivity  has  been  reported 
by  Schwan111  and  Johnson  and  Guy.11  Based  on  the  extrapolation  of  their 
results  for  frequencies  below  1 kHz,  the  realistic  complex  permittivity 
for  a high  water  content  biological  tissue  at  extremely  low  frequencies 
(in  the  range  from  45  Hz  to  76  Hz)  is  considered  to  be  approximately 
e,  = e (10  -jlO  ) farads/m,  where  e„,  the  free  space  permittivity,  is 
8.854  x 10'  farads/m.  This  value  of  is  used  in  the  calculation. 

To  make  the  desired  calculations,  the  following  assumptions  and 
parameters  are  used: 


• The  lattice  spacings  for  the  tissue  culture  tray  of  Figure  3 
are  measured  to  be  9.52,  15.08  and  9.52  mm  along  the  x,  y and 
z directions,  respectively. 

• In  the  tissue  culture  tray,  the  biological  tissue  is  filled 
in  a hollow  cylindrical  hole  whose  geometric  shape  resembles, 
to  a certain  degree,  that  of  a prolate  spheroid.  Hence,  the 
prolate  spheroidal  model  of  biological  tissue  is  used  in  the 
calculation.  For  quantitative  comparison,  calculations  are 
also  made  for  spherical  models  of  biological  tissue. 

• The  biological  tissue  is  embedded  in  a very  thin  layer  of 
plastic  (dielectric  slab)  with  a zero  loss  tangent  and  a 
relative  dielectric  constant  of  2.25  (i.e.,  e9'  = 2.25  e 
and  e2~  = 0). 
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Using  the  prolate  spheroidal  model  of  biological  tissue,  the  induced 
electric  fields  inside  the  tissue  are  calculated  based  on  Equations  (4)-(7). 
The  major  and  minor  radii  of  the  prolate  spheroid,  ra  and  r^,  are  measured 
to  be  6.35  mm  and  3.18  mm,  respectively.  Figure  6 shows  the  calculated 
results  for  a two-dimensional  array.  In  this  case,  the  interaction  con- 
stant, 8,  using  the  expression  of  Equation  (B-ll)  is  calculated  to  be 

3 

-2.25/ira  . The  applied  voltage  between  the  parallel  plates  is  taken  to 
be  1 volt,  and  the  results  are  plotted  as  a function  of  for  two  val- 

ues of  ^2  and  f°r  the  case  °f  = £3’  where  is  the  thickness  of  the 
tissue  culture  tray  and  (or  2,^)  is  the  separation  between  the  tray  and 
the  upper  (or  lower)  plate  of  the  simulator.  Note  that  in  this  figure, 
only  the  magnitude  of  the  field  is  given.  Because  of  the  strong  capacitive 
effect  associated  with  the  tissue  culture  tray  and  the  parallel  plate,  the 
induced  internal  electric  field  E^  and  the  applied  voltage  VQ  are  almost 
90°  out  of  phase. 

Figure  7 illustrates  the  corresponding  results  for  the  spherical  model 

★ 

of  biological  tissue.  The  radius  of  the  sphere,  denoted  by  r , is  chosen 
to  be  about  4 mm  such  that  the  sphere  has  the  same  volume  as  that  of  the 
prolate  spheroid. 

Also  of  interest  to  the  coupling  problem  is  the  three-dimensional  lat- 
tice or  array.  When  two  or  more  two-dimensional  tissue  culture  trays  are 
stacked  together,  the  structure  becomes  a three-dimensional  array.  Hence, 
for  the  same  lattice  spacings  and  object  sizes,  the  three-dimensional  cul- 
ture tray  would  have  a thickness  larger  than  that  of  a two-dimensional 
culture  tray.  Two  values  of  75  mm  and  150  mm,  are  used  in  the  calcula- 
tion. The  results  (using  Equation  (B-7)  of  Appendix  B for  the  interaction 
constant,  8)  are  graphically  presented  in  Figures  8 and  9 for  prolate  spher- 
oidal and  spherical  models  of  biological  tissue,  respectively. 


* 


In  Appendix  A,  the  radius  of  a dielectric  sphere  is  denoted  by  a. 
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Internol  Electric  Field,  |E|j,  V/m 


Fig.  6 INDUCED  ELECTRIC  FIELD  INTERNAL  TO  THE  PROLATE 
SPHEROIDAL  MODEL  OF  BIOLOGICAL  TISSUE  IN  A TWO- 
DIMENSIONAL  TISSUE  CULTURE  TRAY 
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Fig.  8 INDUCED  ELECTRIC  FIELD  INTERNAL  TO  THE 
PROLATE  SPHERICAL  MODEL  OF  BIOLOGICAL 
TISSUE  IN  A THREE-DIMENSIONAL  TISSUE 
CULTURE  TRAY 
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Fig.  9 INDUCED  ELECTRIC  FIELD  INTERNAL  TO  THE 
SPHERICAL  MODEL  OF  BIOLOGICAL  TISSUE 
IN  A THREE  - DIMENSIONAL  TISSUE  CULTURE 
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Based  on  these  data,  we  note  that  the  induced  electric  field  inter- 
nal to  the  tissue  culture  decreases  as  l^/Z^  increases.  The  maximum  field 
strength  occurs  when  the  parallel  plates  are  placed  directly  on  the  tissue 
culture  tray.  As  anticipated,  the  coupling  is  stronger  for  a biological 
tissue  sample  (or  object)  having  a more  elongated  geometric  shape.  VJhen 
compared  to  the  spherical  model,  the  coupling  is  stronger  for  the  prolate 
spheroidal  model  of  the  tissue  sample.  This  is  easily  seen  from  the  re- 
sults presented  in  Figures  6 and  7. 

To  establish  the  validity  of  the  analysis  developed  in  this  study, 
calculations  using  the  circuit  approach  are  also  made.  The  results  are 
compared  with  the  more  exact  solutions,  using  the  interaction  field  ap- 
proach in  Table  1.  Note  that  the  capacitor,  C2,  in  the  equivalent  circuit 
of  Figure  5(a),  in  most  practical  situations,  can  be  neglected  without 
losing  any  accuracy  in  the  result  obtained.  For  example,  the  impedance 
of  (i.e.,  l/jioCg)  is  calculated  to  be  -j  1.62  x 105  ohms  for  the  para- 
meters and  dimensions  of  present  interest.  When  compared  to  resistance 
R2,  which  is  calculated  to  be  about  32.5  ohms,  the  capacitor,  C2  may  be 
ignored.  The  close  agreement  of  the  results  between  the  two  methods, 

TABLE  1 

ELF  ELECTRIC  FIELD  COUPLING  TO  BIOLOGICAL  TISSUE 
UNDER  NMRI  SIMULATOR 


Induced  Internal  Electric  Field  |E.|  (V/m)* 

Rigorous  Interaction 

Field  Approach  for  a 
Two-Dimentional  Periodic 
Array 

Approximate  Circuit 
Approach  Using  the 
Geometric  Model  of 
Figure  5(a) 

Approximate  Circuit 
Approach  Using  the 
Geometric  Model  of 
Figure  5(b) 

1.17  x 10-7 

2.05  x 10-7 

0.71  x lO'7 

★ 

Z = 1 mm,  = 15  mm,  *4  = *-3  = 70  mm 


r = 6.35  mm,  r.  = 3.18  m 

a D 

SR  = 0.349,  Sc  = 0.651 

V = 1 volt. 

0 
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as  indicated  in  the  Tab^e,  suggests  that  the  solutions  developed  in  this 
study  can  be  used  with  confidence  to  calculate  the  induced  electric  field 
inside  the  tissue  culture  used  in  the  NMRI  experiment.  It  also  indicates 
that  the  simple  circuit  approach,  although  admittedly  approximate, provides 
a reasonable  estimate.  In  fact,  the  results  given  in  Table  1 for  the  geo- 
metric models  of  Figures  5(a)  and  5(b)  represent,  respectively,  the  upper 
and  lower  bound  values  for  the  solution  in  question.  The  actual  electric 
field  that  would  be  induced  inside  the  tissue  should  have  a value  smaller 
than  that  predicted  from  the  geometric  model  of  Figure  5(a),  but  larger 
than  that  predicted  from  Figure  5(b). 

Using  the  methods  of  Appendix  A,  it  can  be  shown  that  for  an  external 

field  of  20  V/m,  the  field  internal  to  a biological  body  of  permittivity 

5 8 7 

el  = S(10  -Jl°  ) wou^d  be  approximately  6 x 10  V/m.  To  develop  this 

field  in  tissue  samples  in  a plastic  tray  like  the  one  shown  in  Figure  2, 
requires  5.2  volts  across  the  plates  of  the  simulator  for  the  present 
plate  spacing  (approximately  6 inches).  This  is  for  a single  layer  of 
trays.  For  five  layers  of  trays,  2.1  volts  is  required.  These  calcula- 
tions were  made  using  the  equations  and  graphs  of  Sections  2 and  3.  Linear 
interpolation  may  be  used  to  determine  the  required  voltage  across  the 
plates  for  a number  of  layers  of  trays  between  one  and  five. 

The  tissue  samples  near  the  edges  of  the  trays  will  experience  some- 
what stronger  exposure  than  the  rest  of  the  array.  This  has  not  been  quan- 
tified further. 


I I T RESEARCH  INSTITUTE 


22 


SUMMARY  AND  CONCLUSIONS 


1 


The  Naval  Medical  Research  Institute  (NMRI)  is  performing  an  experi- 
ment to  study  the  possible  effects  of  the  electric  and  magnetic  fields 
associated  with  the  Navy's  proposed  ELF  Communications  System  (Seafarer) 
on  human  tissue  cultures.  This  report  presents  the  results  of  a support- 
ing study  of  the  electric  field  coupling  to  the  tissue  culture  samples. 

Two  analytical  methods  are  developed  in  this  study.  The  first  method, 
which  considers  the  dipole  field  interaction,  is  more  accurate.  The  second 
method  is  very  simple,  using  an  approximate  equivalent  circuit  treatment. 

The  agreement  of  the  results  obtained  from  these  two  methods  is  fairly  good. 
Therefore,  the  resultant  expressions  can  be  used  with  confidence  in  the 
studies  of  ELF  electric  field  coupling  to  biological  tissue  and  in  many 
problem  areas  concerning  the  electromagnetic  biological  effects. 

Observations  and  conclusions  drawn  from  this  study  are: 

(1)  The  ELF  magnetic  field  will  penetrate  the  tissue  cultures 
without  change. 

(2)  As  anticipated,  the  induced  electric  field  internal  to  the 
biological  tissue  in  the  tissue  culture  tray  decreases  as 
the  permittivities  of  the  tissue  and  its  surrounding  medium 
increases. 

(3)  The  ELF  electric  field  internal  to  the  tissue  sample  will  be 
much  less  than  the  applied  field.  For  a single  layer  of  trays, 
the  reduction  factor  is  about  6 x 10'. 

(4)  The  electric  field  coupling  to  the  tissue  also  depends  upon 
the  separation  between  the  tissue  culture  tray  and  the  upper 
and  lower  plates  of  the  simulator.  Maximum  internal  field 
strength  is  achieved  when  the  parallel  plates  are  placed 
right  on  and  beneath  the  culture  tray. 

(5)  Like  the  ELF  electric  field  coupling  to  a single  isolated  bio- 
logical object,  the  coupling  is  stronger  for  objects  having  a 
more  elongated  geometric  shape.  The  coupling  is  strong  when 
the  object  exhibits  a geometric  shape  somewhat  like  a long  and 
thin  circular  cylinder.  On  the  other  hand,  the  coupling  is 
small  for  geometric  shape  somewhat  like  a flat  circular  disk. 
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(6)  From  the  limited  numerical  results  obtained  in  this  study  (see 
Appendix  C),  it  appears  that  the  array  spacing  and  the  object 
size  also  affect  the  electric  field  coupling.  However,  the 
effect,  which  primarily  arises  from  the  variation  in  the 
effective  permittivity  of  the  tissue  culture  tray  for  differ- 
ent array  spacing  and  object  size,  is  not  significant  at 
extremely  low  frequencies. 

t (7)  Although  sophisticated  field  theory  has  been  required  to 

analyze  the  complex  geometry  used  in  this  experiment,  an 
adequate  estimate  of  the  field  coupling  can  be  obtained 
using  a simple  circuit  model.  This  model  also  aids  in  the 
understanding  of  the  coupling  mechanism. 
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1.  INTRODUCTION 


The  electric  field  coupling  to  and  penetration  of  dielectric  objects 
depends  on  the  frequency  of  the  electric  field,  the  complex  permittivity 
(i.e.,  the  dielectric  constant  and  conductivity)  of  the  object  and  its 
surrounding  medium,  and  the  size  and  shape  of  the  object.  At  extremely 
low  frequencies  (ELF)  the  simplifying  assumption  is  made  that  the  object 
is  small  compared  to  wavelength  and  then  only  the  shape  of  the  object  and 
the  permittivities  involved  are  parameters.  Of  interest  are  the  electric 
field  internal  to  the  object,  , the  polarization,  P,  and  the  polariza- 
bility, ag,  of  the  object. 

The  simple  geometric  shapes  to  be  considered  here  are  illustrated  in 
Figure  A-l.  The  cylinders  of  Figures  A-l(b)  and  A-l(c)  are  assumed  to  be 
very  thin  and  the  slab  of  Figure  A-l(d)  is  also  assumed  to  be  very  thin. 

The  prolate  (or  oblate)  spheroid  is  formed  by  the  rotation  of  an  ellipse 
about  its  major  (or  minor)  axis.  In  each  case  the  permittivity  of  the 
outside  medium  is  designated  ^ and  that  of  the  inside  medium  is  e^.  These  per- 
mittivities are  complex.  For  example, 

E1  = (A~1} 

where 

= permittivity 


Oj  = conductivity 
to  = 2nf 
f = frequency. 

2.  INTERNAL  ELECTRIC  FIELD 

The  induced  field  within  the  dielectric  object  when  exposed  to  an  ELF 
electric  field  may  be  obtained  for  the  simple  geometric  shapes  of  Figure  A-l. 
The  results  are  presented  in  Table  A-l. 
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(a ) Sphere 


(b)  Vertical  Cylinder 
h»  2a 
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(c  ) Horizontal  Cylinder 
h»2a 


Fig.  A- 1 DIELECTRIC  OBJECTS  OF  SIMPLE  GEOMETRIC 
SHAPE 
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The  static  solution  for  a dielectric  sphere  has  been  obtained  by 

Stratton  (pp  205-206).  The  sphere  has  been  used  by  many  investigators 

12  13 

as  an  approximate  geometric  model  for  human  baby.  ’ Note  that  the 

internal  field  is  polarized  along  the  same  direction  as  the  electric  field 
impressed  on  the  object.  The  internal  field  is  uniform  and  is  independent 
of  sphere  size. 

For  the  long  thin  cylinder  of  Figure  A-l(b)  with  impressed  electric 
field  along  its  axis,  the  continuity  of  the  tangential  electric  field  along 
the  surface  boundary  requires  that  ~ EQ.  Again,  this  approximation  is 
independent  of  the  cross-sectional  shape  of  the  cylinder  and  is  valid  when 
the  cylinder  is  long  and  thin. 

In  Figure  A-l(c),  the  external  electric  field  is  normal  to  the  axis 
of  the  thin  cylinder.  The  result  given  in  Table  A-l  is  the  static  solution 

14  1 

obtained  by  Smyth. 

For  the  thin  slab  with  the  electric  field  normal  to  it,  the  continuity 
of  the  normal  component  of  electric  flux  requires  that  C2^0  = ei^i’ 

The  prolate  spheroid  shown  in  Figure  A-l(e)  is  a symmetric  ellipsoid 
which  is  generated  by  revolving  an  ellipse  about  its  major  axis.  The  oblate 
spheroid,  as  shown  in  Figure  A-l(f),  is  formed  by  the  revolution  of  an 
ellipse  about  its  minor  axis.  A very  general  static  solution  for  the  case 

g 

of  an  ellipsoid  has  been  developed  by  Stratton  (pp  211-213),  but  the  re- 
sultant expression  is  in  terms  of  elliptical  integrals  which  usually  require 
approximate  numerical  treatment.  A conclusion  of  great  practical  importance 
that  can  be  reached  from  Stratton's  results  is  that  if  the  applied  field 
is  initially  uniform  and  parallel  to  one  of  the  ellipsoid  axes,  then  the 
resultant  internal  field  within  the  ellipsoid  is  also  uniform  and  parallel 
to  that  same  axis.  In  Table  A-l,  we  list  the  solutions  for  the  prolate  spheroid 
and  the  oblate  spheroid.  These  low  frequency  results  are  obtained  for  these 
two  non-spherical  spheroids  using  the  method  of  separation  of  variables.^ 

The  parameters  u and  IT  for  the  prolate  spheroid  and  the  oblate  spheroid  are 
given  by 


i 
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u 


21-1/2 


• ■ Hji  1 ^ 


21-1/2 


(A-2) 


where  2a  and  2b  are  the  major  and  minor  axes  of  the  spheroid,  respectively. 
Since  the  limit  as  b-*-a  corresponds  to  the  limit  of  u-«°  (or  it  can  be 
shown  that  under  this  limiting  case,  the  induced  internal  electric  field 
becomes  E.  = 3E0/(c^/e2+2)  for  both  prolate  and  oblate  spheroids.  As  ex- 
pected, this  is  identical  to  the  solution  obtained  for  a dielectric  spheroid. 

The  field  boundary  conditions  require  the  magnitude  of  the  tangential 
internal  electric  field  to  approach  that  of  the  tangential  external  field. 

The  field  boundary  conditions  also  require  the  normal  component  of  the  in- 
ternal field  to  be  that  of  external  fields  reduced  by  a factor  of  z^/^i 
which,  in  general,  is  a small  quantity.  Hence,  as  can  be  seen  from  Table  A-l, 
the  coupling  of  ELF  electric  field  into  the  object  is  the  strongest  for  the 
very  thin  cylinder  as  shown  in  Figure  A-l(b),  because  the  incident  electric 
field  is  tangential  to  most  of  its  surface  boundary.  For  the  prolate  spher- 
oid of  Figure  A-l(e),  the  incident  electric  field  is  tangential  to  the  long 
major  axis  of  the  spheroid.  Thus,  the  electric  field  coupled  into  the  pro- 
late spheroid  is  the  next  strongest.  When  compared  to  prolate  spheroid,  the 
oblate  spheroid  provides  a weaker  coupling,  because  for  the  same  b/a  ratio, 
the  prolate  spheroid  represents  a more  elongated  geometry  than  the  oblate 
spheroid.  For  the  symmetric  dielectric  sphere,  the  internal  induced  field 
strength  is  smaller  than  that  of  a prolate  spheroid,  but  is  larger  than  the 
field  coupled  into  an  oblate  spheroid.  The  geometry  of  a horizontal  circular 
cylinder  can  provide  a relatively  small  electric  field  coupling,  but  the 
coupling  is  the  smallest  for  the  thin  slab  of  Figure  A-l(d),  because  the 
external  impressed  field  is  perpendicular  to  most  of  its  surface  boundary. 

This  demonstrates  that  the  ELF  electric  field  coupling  to  dielectric  objects 
is  strongly  dependent  on  the  geometric  shape  of  the  object,  but  not  dependent 
on  the  size  of  the  object  which  is  much  smaller  than  wavelength  at  extremely 
low  frequencies. 
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3.  POLARIZATION  AND  POLARIZABILITY 

CONSTANT  OF  DIELECTRIC  OBJECTS 

A characteristic  of  all  dielectric  materials  is  their  ability  to  store 
electric  energy  under  the  influence  of  an  electric  field.  This  storage  of 
energy  takes  place  by  a shift  in  the  relative  positions  of  the  internal  posi- 
tive and  negative  charges  against  normal  molecular  and  atomic  forces.  Under 
the  influence  of  a given  electric  field,  the  amount  of  the  charge  displace- 
ment differs  for  different  dielectric  materials.  Electrical  dipole  or  simply 
a dipole  is  the  name  given  to  two  point  charges  of  equal  magnitude  and 
opposite  sign  separated  by  a small  distance.  The  term  "polarization"  or 
(dipole  moment)  is  used  to  characterize  the  nature  of  dielectric  material 
and  its  geometric  shape  when  exposed  to  an  external  electric  field. 

Consider  a dielectric  object  embedded  in  a homogeneous  and  isotropic 
medium.  The  permittivities  of  the  dielectric  object  and  its  surrounding 
medium  are  designated  by  and  respectively.  Under  the  influence  of  a 
uniform  external  field,  EQ,  the  polarization  of  the  object,  p,  defined  by 
Stratton8  (pp  183  and  206)  is  given  by 

P = (ei“e2 ) Ei 

As  illustrated  in  Table  1,  the  internal  induced  field,  E- , is  related  to  the 
external  field  EQ.  Hence,  the  polarization  of  a dielectric  object  can  be 
expressed  in  terms  of  the  external  field,  EQ.  As  an  example,  the  polariza- 
tion of  a dielectric  sphere,  from  Equation  (A-3),  is 


P 


3e, 


£r£2 

+ 2E;  fco 


'1 


(A-4) 


Since  the  polarization  is  also  defined  as  the  dipole  moment,  p,  per  unit 
volume.  Equation  (A-4)  gives 


P 


4na2  e. 


I ej-e2 
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(A-5) 
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The  relationship  between  polarization  P (or  dipole  moment  p)  and  the 
impressed  external  field,  EQ,  is,  of  course,  a function  of  the  type  of  di- 
electric material.  For  isotropic  material  which  is  the  case  of  present 
interest,  P (or  p)  and  EQ  are  linear  related  and  are  always  parallel 
regardless  of  the  orientation  of  the  external  field,  The  constant  which 
characterizes  this  relationship  between  p and  EQ  is  termed  the  polarizabil- 
ity, ag,  of  the  dielectric  object. 

P = aes2EQ  (A-6) 


To  illustrate  this,  consider  the  dielectric  sphere  of  Figure  A-l(a). 
A comparison  between  Equations  (A-5)  and  (A-6)  gives  the  polarizability 
constant  for  the  dielectric  sphere 


a 


e 


4iTa3 


where  £ = 


(A-7) 


Figure  A-2  illustrates  the  polarizability  for  a dielectric  sphere  with 

a real  permittivity  (i.e.,  e,"  = 0).  In  the  figure,  the  polarizability, 

3 * 

ag,  is  normalized  to  4ira  , where  a is  the  radius  of  the  sphere.  It  can  be 
seen  that  a becomes  zero  if  £ is  taken  to  be  unity.  This  is  expected  be- 
cause  when  no  dielectric  object  exists  in  the  medium.  It  is  also 

of  interest  to  consider  the  case  of  perfect  conducting  objects.  Noting  that 
the  ductivity  of  a dielectric  sphere  with  a complex  permittivity 
is  given  by  , a dielectric  object  with  z^"  approaching  <»  would  be 

equivalent  to  a perfect  conducting  body.  The  polarizability  of  Equation  (A-7) 

3 

under  this  limiting  case  of  z^" -*■<*>  gives  = 4rra  , the  polarizability  of 
a conducting  sphere.  This  result  agrees  with  those  used  in  the  theory  of 
artificial  dielectric.6,7 


Using  Equations  (A-3)  and  (A-6)  the  polarization  P and  the  polarizabil- 
ity ag  for  the  dielectric  objects  of  Figure  A-l  can  be  obtained.  The  resul- 
tant expressions  are  summarized  in  Table  A-2.  These  simple  expressions  are 
useful  in  studies  of  artificial  dielectric,  waveguide  array  and  the  inter- 
action field  of  dielectric  objects  in  a periodic  lattice. 
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Normalized  Polarizability,  a 


APPENDIX  B 

STATIC  INTERACTION  FIELDS  IN  A PERIODIC  LATTICE 


I 
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1.  INTRODUCTION 


The  electric  field  coupling  to  dielectric  objects  can  be  discussed  in 
terms  of  polarization  (or  dipole  moment),  because,  under  the  influence  of 
an  external  field,  charges  are  induced  in  the  object  and  are  displaced  so 
as  to  set  up  an  induced  field  that  will  eventually  satisfy  the  boundary 
conditions  on  the  dielectric  surface.  Each  dielectric  object  generates  a 
dipole  field  and,  thus,  exhibits  a dipole  moment.  When  a great  number  of 
dielectric  objects  are  embedded  in  a homogeneous  and  isotropic  medium, 
interaction  occurs  between  all  of  the  dipole  fields  produced  by  the  objects. 
In  this  Appendix,  a suitable  analytical  method  is  presented  to  determine 
the  interaction  fields  for  the  dielectric  objects  in  a periodic  lattice. 

Of  present  interest  is  the  interaction  of  electric  fields  at  extremely  low 
frequencies  (ELF).  However,  the  results  can  be  easily  extended  to  include 
the  magnetic  field  coupling  and  the  case  of  plane  electromagnetic  fields. 

2.  INTERACTION  FIELDS  FOR  DIELECTRIC  OBJECTS  IN  A 

THREE-DIMENSIONAL  PERIODIC  LATTICE 

The  periodic  lattice  considered  here  is  shown  in  Figure  B-l,  where  the 
spacings  between  objects  are  designated  by  a,  b,  and  c along  the  x,  y,  and 
z axes,  respectively.  At  extremely  low  frequencies,  only  the  static  inter- 
action field  is  of  interest,  because  the  size  of  the  object  and  the  spacings 
between  the  objects  are  sufficiently  small  compared  to  wavelength. 

To  obtain  the  static  interaction  field  for  the  lattice  of  Figure  B-l, 
each  dielectric  object  may  be  represented  by  its  equivalent  dipole  moment,  p. 
The  dipole  moment  (and  the  electric  field  coupling)  for  dielectric  objects 
of  simple  geometric  shapes  has  been  discussed  in  Appendix  A.  For  objects 
that  are  symmetrically  distributed,  the  image  theory  of  perfect  conducting 
planes  may  be  employed.  It  can  be  seen  from  Figure  B-l  that  conducting 
planes  may  be  inserted  into  the  lattice  structure  at  y = _+  b/2  without  dis- 
turbing the  field  distribution.  Hence,  by  using  image  theory,  the  potential 
distribution  of  the  lattice  of  Figure  B-l  may  be  obtained  by  first  consid- 
ering the  potential  for  a y-directed  dipole  located  at  y = 0 between  the 
conducting  plates  at  y = + b/2  as  shown  in  Figure  B-2,  and  then  summing  the 
potential  distributions  for  all  the  dipoles  placed  at  x = na,  y = 0 and 
z = sc  for  n = 0,  + 1,  + 2 • ••  and  s = 0,  + 1,  + 2"*. 
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Fig.B-l.  DIELECTRIC  OBJECTS  IN  A THREE-DIMENSIONAL 
PERIODIC  LATTICE 
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RIBUTION  FOR  AN  ELECTRIC 
N TWO  CONDUCTING  PLATES 
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The  dipole  at  the  origin  may  be  considered  as  a positive  charge 
q 6(x)  6(z)  6(y-A)  atx  = z = 0,  y = A and  a negative  charge  -q  6(x) 

6(z)  <5(y+A)  at  x = z = 0,  y = -A,  where  6(x-xQ)  represents  the  delta 
function  which  is  defined  to  be  equal  to  zero  at  all  points  of  x except 
at  x ; xQ  where  it  becomes  infinite.  The  potential  function  for  the 
lattice  of  Figure  B-2  is  governed  by  the  static  Poisson's  equation  which 
in  cylindrical  coordinates  (p,0,y)  is  given  by 


A , 1 JL 
2 

9y  p 3p 


pi 

= -^-6(x)6(z) 

6(y-A)-6(y+A) 

l 9p 

e2 

■ 

(B-l) 
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The  solution  of  Equation  (B-l)  has  been  obtained  by  Collin  and  the 
resultant  expression  is  given  by  Equation  (B-3)  below.  The  mathematical 
treatment  is  briefly  presented  as  follows. 


The  potential  function  $ is  independent  of  6.  It  can  also  be  seen  from 

Figure  B-2  that  is  an  add  function  of  y and  must  vanish  on  the  conducting 

planes  at  y = + b/2.  The  solutions  of  Equation  (B-l)  having  a variation  of 

sin(2(mry/b)  with  respect  to  y,  are  the  functions  of  the  modified  Bessel  func 

tions  of  the  first  and  second  kind.  Only  the  second  kind  of  the  modified 

Bessel  function  K (2imrp/b)  i s allowed,  because  it  satisfies  the  radiation 
0 22  1/2 

condition  that  $>  -*■  0 as  p = (x  +z  ) ' becomes  large.  Hence,  the  general 
solution  for  the  potential  function  0 may  be  expressed  as 


00 


(B-2) 


where  the  constant  Am  is  determined  from  the  dipole  source  as  p -*■  0.  Noting 

that  the  dipole  moment  p is  defined  as  p = 2qA  as  A approaches  zero,  the 

2 

constant  Am  is  2mp/e2b  , where  is  the  permittivity  of  the  medium  in  which 
dielectric  objects  are  embedded.  Summing  up  the  potential  distributions  due 
to  the  dipoles  along  the  x and  z axes  and  noting  that  the  electric  field  is 
related  to  the  potential  function  by  E = -84>/9y,  the  resultant  expression 
for  the  interaction  field  due  to  all  dipoles  including  the  one  at  the  origin 
is  given  by 
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(B-3) 


where  the  prime  in  the  summation  over  m indicates  the  omission  of  the  term 
m = n = s = 0. 


The  modified  Bessel  function  of  the  second  kind,  KQ(x),  has  the  property 
that  as  x increases,  K0(x)  decreases  rapidly.  This  is  evidenced  from  Fig- 
ure B-3  in  which  K (x)  is  plotted  as  a function  of  x.  Using  the  fact  that 

o 16 

Kq(x)  decays  exponentially  at  a rapid  rate  and  that 


(B-4) 


the  interaction  field  at  the  center  of  the  dielectric  object  (at  y=0),  from 
Equation  (B-3),  may  be  approximated  by  the  following  simple  form 


1.201 
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(B-5) 


In  most  cases,  Equation  (B-5)  provides  sufficient  accuracy  for  the  s^ution, 
especially  when  a and  c do  not  differ  greatly  from  b.  A simple  numerical 
computation  indicates  that,  if  a = b = c,  the  calculated  value  for  the  largest 
term  neglected  in  Equation  (B-5)  is  less  than  0.2  percent  that  of  the  second 
term  of  Equation  (B-5).  The  expression  of  Equation  (B-5)  represents  the 
interaction  field  at  the  center  of  the  dielectric  object  arising  from  the 
interaction  of  all  the  dipole  fields  generated  in  the  periodic  lattice. 


The  relationship  between  the  interaction  field,  Ey^ , and  the  dipole 
moment  of  a single  (isolated)  dielectric  object  is  linear,  as  is  clearly 
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indicated  from  Equation  (B-5).  The  constant  which  characterizes  this  rela- 


tionship is  the  interaction  constant.  The  interaction  constant,  defined 

6,9 


as 


6 = E 


y1 


(B-6) 


has  been  extensively  used  in  the  theory  of  artificial  dielectrics.  From 
Equations  (B-5)  and  (B-6),  the  interaction  constant,  0,  for  the  three-dimen- 
sional periodic  lattice  of  Figure  1 is  given  by 
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(B-7) 


3. 


INTERACTION  FIELDS  FOR  DIELECTRIC  OBJECTS 
IN  A TWO-DIMENSIONAL  PERIODIC  LATTICE 


The  two-dimensional  periodic  lattice  considered  here  is  illustrated 
in  Figure  B-4,  where  every  dielectric  object  in  the  lattice  has  been  re- 
placed by  its  equivalent  dipole  moment  p.  To  obtain  the  static  interaction 
field  of  this  lattice,  we  note  that  the  three-dimensional  interaction  field 


of  Equation  (B-3)  may  be  expressed  in  an  alternative  form  as 
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2 (mb ) 2 - (na )2  - (sc)2 


m=  —co  n=  -°°  s=  -oo 

In  obtaining  Equation  (B-8),  the  following  Poisson  summation  formula0’1' 


(mb)2  + (na)2  + (sc)2 
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has  been  appropriately  applied  where  F(u)  is  the  Fourier  transform  of  f(t). 
A comparison  between  Figure  B-4  and  Figure  B-l  reveals  that  the  inter- 


action field  for  the  two-dimensional  lattice  may  be  obtained  from  Equation 
(B-8)  by  omitting  the  summation  over  m.  Again,  using  the  Poisson  summation 


formula  and  the  approximation  that  the  modified  Bessel  function  Kq(x)  de- 


creases rapidly  as  x increases,  it  can  be  shown  that  the  interaction  field 
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J. 


B-4  DIPOLE  REPRESENTATION  FOR  DIELECTRIC  OBJECTS 
IN  A TWO-DIMENSIONAL  LATTICE 


at  the  center  of  the  dielectric  object  for  the  two-dimensional  lattice  of 
Figure  B-4  is  given  by 


APPENDIX  C 

DETERMINATION  OF  THE  FIELDS  INDUCED  INSIDE  THE  BIOLOGICAL  OBJECT 
FOR  PARALLEL-PLATE  SIMULATOR  EXPERIMENTATION 
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1.  INTRODUCTION 

An  analytical  approach  is  described  in  this  Appendix  to  the  problem 
of  determining  the  field  induced  inside  the  biological  object  (tissue)  for 
the  parallel-plate  simulator  experimentation  at  NMRI.  The  presentation  is 
based  on  the  analysis  being  generally  applicable  and  not  dependent  on  the 
specific  details  of  the  simulator  system.  The  result  can  be  used  to  cal- 
culate the  electric  field  internal  to  the  biological  object  for  arbitrary 
dimensions  of  the  simulator  and  for  biological  objects  of  any  geometric 
shape,  dimension  and  dielectric  property.  With  the  polarizahi 1 ity  constant 
of  the  dielectric  object  (see  Appendix  A)  and  the  static  interaction  constant 
in  the  periodic  lattice  (see  Appendix  B)  obtained,  this  analytical  develop- 
ment is  accomplished  by  using  the  theory  of  artificial  dielectric.  Specific- 
ally, we  wish  to  calculate  the  internal  field  within  the  biological  objects 
in  a periodic  lattice  in  terns  of  the  voltage  applied  between  the  two  plates 
of  the  simulator  as  well  as  the  physical  characteristics  and  the  electrical 
parameters  of  the  biological  object. 

In  this  discussion,  a realistic  model  for  the  simulator  which  neglects 
the  distortion  of  the  fringing  field  between  the  parallel  plates  has  been 
adopted  such  that  the  field  external  to  the  periodic  lattice  is  always  uni- 
form. This  condition  implies  that  the  physical  dimension  of  the  parallel- 
plate  simulator  is  larger  than  that  of  the  periodic  lattice. 

2.  EFFECTIVE  PERMITTIVITY  OF  A PERIODIC  LATTICE 

2 . 1 Theoretical  Basis 

The  geometrical  construction  for  the  paral lel-plate  simulator  consid- 
ered here  is  illustrated  in  Figure  C-l.  The  biological  objects  under  test 
are  embedded  in  a homogeneous  and  isotropic  dielectric  slab  in  the  form  of 
a tissue  culture  tray.  The  permittivities  of  the  biological  objects  (tissue) 
and  its  surrounding  medium  are  designated  by  and  respectively.  The 
periodic  lattice  is  exposed  to  a uniform  electric  field  of  very  low  frequency 
which,  in  the  case  of  present  interest,  is  generated  by  a paral lel-plate. 

The  periodic  lattice  with  a height  of  H,>  is  located  between  the  parallel 
plates  at  a distance  from  the  upper  plate  and  from  the  lower  plate. 


IIT  RESEARCH  INSTITUTE 

46 


Porallel  Plates  L 


k. 


*1 


i v- 


Periodic  Lattice 


1 


y 


Fig.C.I  A PERIODIC  LATTICE  EXPOSED  TO  A UNIFORM  ELECTRIC  FIELD 
GENERATED  BY  A PARALLEL-PLATE  SIMULATOR 


Under  the  influence  of  EQ,  electric  charges  are  induced  in  the  bio- 
logical objects,  Each  object  exhibits  a dipole  field.  As  a result,  inter- 
action occurs  in  the  periodic  lattice  between  all  the  dipole  fields  generated 

by  the  biological  objects,  The  net  combined  effect,  as  has  been  used  in  the 

6 7 

study  of  an  artificial  dielectric,  ’ may  be  described,  in  terms  of  effec- 
tive permittivity  ee>  by 


D * ea  En 
e p 


(C-l) 


where  Ep  denotes  the  electric  field  that  would  be  induced  in  the  periodic 
lattice. 

Based  on  the  theory  of  artificial  dielectric,  the  effective  permittivity 
for  the  periodic  lattice  of  Figure  C-l  may  be  expressed  as& 


\ Nae 

= e,  <1  + , — — q 

2 ) l-a„ 6 


(C-2) 


where  ag  is  the  polarizability  constant  which  characterizes  the  nature  of 
the  polarization  under  impressed  electric  fields  in  biological  objects,  N 
is  the  number  of  objects  per  unit  volume  and  B is  the  interaction  constant. 
In  Appendix  A,  the  polarizability  constant  has  been  discussed  in  detail  for 
dielectric  objects  of  various  geometric  shapes.  The  interaction  constant  B 
of  a three-dimensional  (or  two-dimensional)  periodic  lattice  has  also  been 
presented  in  Appendix  B.  The  value  of  N in  Equation  (C-2)  is  equal  to 
(abc)  * for  the  case  of  Figure  C-l,  where  a,  b,  and  c denote  the  spacings 
between  the  biological  objects  along  the  directions  of  x,  y and  z axes, 
respectively. 


Provided  eg  is  known,  the  electric  field  internal  to  the  periodic 
lattice  may  then  be  obtained  from 


£oEo 


(C-3) 


where 
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(C-4) 
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and  Vo  is  the  applied  voltage  between  the  plates,  In  Equations  (C-3)  and 
(C-4)  we  have  used  the  relation  VQ  = U^+£3)Eo+2.gEp  anc*  ^e  conservation  of 
electric  flux  density  (i.e.,  the  normal  component  of  D must  be  continuous 
at  the  surface  boundary  of  the  periodic  lattice). 

To  illustrate  the  physical  insight,  it  appears  constructive  to 
consider  some  special  cases.  For  example,  if  the  biological  object  consid- 
ered here  has  the  same  permittivity  as  that  of  its  surrounding  medium 
= £g,  then  the  periodic  lattice  becomes  a dielectric  slab,  as  shown  in 
Figure  C-2(a).  Under  this  condition,  £ = e^/eg  = 1 anc*  ae  = 0 (see  Appen- 
dix A),  we  then  have  from  Equation  (C-2)  that  eg  = eg  which,  in  turn,  gives 
^p’a  = eoEo^e2"  as  expected,  is  the  solution  of  a dielectric  slab 

with  a permittivity  Eg  immersed  in  a uniform  impressed  electric  field,  EQ. 
Likewise,  the  electric  field  internal  to  the  periodic  lattice  becomes 
^p’b  = £oEo/el  ^ t^ie  tn°l°9ica1  objects  with  a permittivity  e^,  shown  in 
Figure  C-2(b),  are  embedded  in  a dielectric  having  the  same  permittivity. 
Evidently,  Ep,a  is  larger  than  Ep>b>  if  Therefore,  we  would  ex- 

pect that  £0  would  have  a value  somewhat  larger  than  eg  but  smaller  than 
e,  and  that 

1,  p,b  p p,a 

2.2  Illustrative  Example 

As  can  be  seen  from  Equation  (C-2),  the  effective  permittivity 
of  a periodic  lattice  depends  upon  many  parameters  which  include  the  permit- 
tivities of  the  biological  object  (ej)  and  its  surrounding  medium  (eg),  the 
size  and  shape  of  the  object  and  the  interaction  constant.  As  an  illustra- 
tive example,  let  us  consider  a three-dimensional  periodic  lattice  of  Fig- 
ure C-l.  The  lattice  consists  of  a large  number  of  biological  objects  which 
are  embedded  in  a dielectric  (or  plastic)  slab.  The  geometry  of  the  object 
is  taken  to  be  a prolate  spheroid.  In  this  example,  we  use 

a = 3/8  inch,  the  lattice  element  spacing  along  the 
direction  of  x-axis; 

b = 19/32  inch,  the  lattice  element  spacing  along  the 

direction  of  y-axis; 

c = 3/8  inch,  the  lattice  element  spacing  along  the 
direction  of  z-axis; 


j 
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1/4  inch,  the  major  radius  of  the  prolate  biological  object; 
1/8  inch,  the  minor  radius  of  the  prolate  biological  object. 


I 


When  exposed  to  a uniform  ELF  electric  field,  the  interaction  constant, 

3 

from  Equation  (B-7)  in  Appendix  B,  is  obtained  to  be  8 a -0 . 62/ irb  . The 
insertion  of  B in  Equation  (C-2)  results  in  the  effective  permittivity  for 
the  periodic  lattice 


e2 


1 + 0.211 
1 + 0.015 


(C-5) 


where 


A = jl+(?-l)(uI 2-l)[u  coth-1  (u)-l]j 


(C-6) 


u = 


r ?i 

C\J 

r— 1 

P^tAa)  J 

Figure  C-3  shows  the  normalized  effective  permittivity  (with  respect 
to  e 2 ) obtained  from  Equations  (C-5)  and  (C-6).  The  normalized  effective 
permittivity  is  plotted  as  a function  of  the  ratio  £ = 

If  > £2,  the  effective  permittivity  of  the  periodic  lattice,  as 
mentioned  previously,  always  has  a value  larger  than  £2  but  smaller  than 
£j.  This  is  evidenced  from  Figure  C-3.  The  effective  permittivity  for  the 
periodic  lattice  at  C = 20,  for  example,  has  a value  which  is  about  1.99 
times  larger  than  ep,  but  is  approximately  10  times  smaller  than  e^.  When 
Ej  = £2,  no  biological  object  exists  in  the  lattice.  In  this  case.  Equa- 
tion (C-4)  as  well  as  Figure  C-3,  as  expected,  predict  that  ce  = = £2* 


3.  INDUCED  ELECTRIC  FIELD  INTERNAL  TO  THE 

BIOLOGICAL  OBJECTS  IN  A PERIODIC  LATTICE 

3 . 1 General  Formulation  of  the  Problem 

When  the  periodic  lattice  is  exposed  to  a uniform  electric  field  EQ, 
electric  charges  are  induced  in  the  biological  objects.  The  charges  in  the 
objects  are  displaced  to  set  up  an  induced  field  that  will  eventually  satisfy 
the  boundary  conditions  at  the  object  surface.  Each  object  thus  exhibits  a 
dipole  field  and,  hence,  interaction  occurs.  At  extremely  low  frequencies 
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Fig.  C-  3 VARIATION  IN  EFFECTIVE  PERMITTIVITY 
FORA  THREE-DIMENSIONAL  PERIODIC 
LATTICE 
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where  the  wavelength  is  long  compared  to  the  dimension  of  the  object,  the 
electric  field  acting  to  polarize  any  object,  therefore,  may  be  assumed  to 
be  a uniform  field  equal  to  the  field  at  the  center  of  the  object.  The 
field,  from  the  considerations  of  the  symmetric  geometry  of  Figure  C-l,  has 
only  a y-directed  component.  The  total  combined  effect  of  all  the  biological 
objects  in  the  periodic  lattice  is  to  produce  a new  average  dipole  polariza- 
tion, P,  for  the  medium.  This  basic  mechanism  may  be  best  described  by  the 
Lorentz  theory  in  which  the  electric  flux  density,  D,  (given  by  Equation 
(C-l))  consists  of  two  parts;  the  first  part  relates  the  electric  field  to 
the  medium  permittivity,  e2,  and  the  second  part  is  called  the  dipole  polar- 
ization. That  is 


£2Ep  + P 


(C-7 ) 


From  this  expression  as  well  as  Equations  (C-l)  and  (C-2),  it  becomes 
straightforward  to  obtain  the  new  average  dipole  moment,  p,  (i.e.,  P/N  where 
N = abc) 


l-age 


(C-8) 


Inserting  Ep,  given  by  Equation  (C-3),  into  Equation  (C-8),  the  dipole  moment, 
p,  may  also  be  expressed  in  terms  of  VQ,  the  applied  voltage  between  the  upper 
and  lower  plates  of  the  simulator,  as 


1-a  6 


e6/  ([^2+^i+!l3^ee^eo^] 


(C-9) 


To  determine  the  electric  field  strength  induced  internal  to  the  bio- 
logical objects  in  a periodic  lattice,  we  note  that  the  dipole  moment,  p, 
as  defined  by  Stratton,  is  related  by  the  internal  electric  field,  E^, 
by 


M£rE2)Ei 


(C-10) 


where  ky  is  the  volume  of  the  object,  since  the  polarization  is  defined  as 
the  dipole  moment,  p,  per  unit  volume.  Using  the  expression  of  Equation  (C-9) 
which  is  the  total  dipole  moment  including  the  effect  of  (dipole)  field 
interactions,  we  have 


I I T RESEARCH  INSTITUTE 


(C-ll) 


Ei  = 


[1-a  6) 


V0/(5-D 


e6)kv)  c[£2+Ui+£3)(ee/E0)] 


for  the  electric  field  induced  inside  the  biological  object  when  the  periodic 
lattice  is  exposed  to  a uniform  electric  field  generated  by  a parallel  plate 
voltage,  VQ. 

From  the  polarizability  constants  listed  in  Table  A-2  in  Appendix  A, 
the  expression  of  Equation  (C-ll)  reveals  that  if  Cj  = z^  (i.e.  £ = 1, 
ee  = e1  = z^)-,  the  induced  electric  field  internal  to  the  biological  object 
for  any  geometric  shape  becomes  E..  = v0/[S'2+^i+*'3^el^eo]’  which>  as  exPec~ 
ted,  is  the  result  obtained  from  a dielectric  slab  situated  between  the 
parallel  plate. 

3.2  Parametric  Results 

In  this  section,  sample  calculations  for  the  induced  electric  field 
internal  to  the  biological  objects  in  a periodic  lattice  are  presented.  The 
calculation  is  based  on  Equation  (C-ll).  Note  that  the  value  of  , the 
permittivity  of  the  biological  object,  employed  here  is  for  exploratory 
parametric  calculations,  and  is  different  from  that  used  in  the  text  where 
the  permittivity  for  a high  water  content  biological  tissue  is  used. 

We  first  consider  the  three-dimensional  lattice.  Calculations  have 
been  carried  out  for  dielectric  objects  of  three  geometric  shapes;  prolate 
spheroid,  sphere  and  finite  circular  cylinder.  In  each  case,  the  dimension 
of  the  dielectric  object  is  properly  chosen  such  as  to  have  the  same  volume. 
Note  that  the  resultant  expression  for  the  polarization  constant  obtained  in 
Appendix  A for  a long  thin  cylinder  is  used  for  the  finite  cylindrical  di- 
electric object.  The  results  for  the  finite  cylinder  model  of  dielectric 
objects,  therefore,  are  approximate,  but  are  included  here  for  comparison. 

In  this  calculation,  the  following  notations  and  parameters  are  used: 

a = 3/8  inch,  the  lattice  element  spacing  along 
the  direction  of  x-axis$ 

b = 19/32  inch,  the  lattice  element  spacing  along 

the  direction  of  y-axis; 

c = 3/8  inch,  the  lattice  element  spacing  along 
the  direction  of  z-axis-, 
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1/4  inch,  the  major  radius  of  the  prolate  spheroidal 
dielectric  object; 


rb  = 1/8  inch,  the  minor  radius  of  the  prolate  spheroidal 

dielectric  object; 

rQ  = 5/32  inch,  the  radius  of  the  dielectric  sphere. 

r = 1/8  inch,  the  radius  of  the  finite  cylindrical 
dielectric  object; 

h = 7/16  inch,  the  height  of  the  finite  cylindrical 
dielectric  object; 

JL  = spacing  between  the  periodic  lattice  and  the  upper  plate 
of  the  simulator; 

JL,  = height  of  the  periodic  lattice; 

JU  = spacing  between  the  periodic  lattice  and  the  lower  plate 
of  the  simulator; 

L = separation  between  the  upper  and  lower  plates  of  the 
simulator; 

e,  = 81  e , the  permittivity  of  the  dielectric  object,  where 

l 0 n 

eQ  = 10  /36u  farads/m; 

= 2.25  eQ  farads/m,  the  permittivity  of  the  dielectric  slab 
in  which  dielectric  objects  are  embedded; 

£ - ej/£2’ 

VQ  = voltage  applied  between  the  plates  of  the  simulator. 

The  calculated  results,  obtained  for  the  case  of  = L/4,  5^2  = L/2 , 

are  graphically  presented  in  Figures  C-4  to  C-6  for  the  sphere,  prolate  spher- 
oid and  finite  cylinder  dielectric  objects,  respectively.  In  the  figure,  the 
induced  internal  electric  field  is  normalized  to  (VQ/L). 

We  now  consider  the  two-dimensional  periodic  lattice.  The  dielectric 
objects  are  embedded  in  a dielectric  slab  in  the  x-z  plane  (see  Figure  C-l) 
with  lattice  spacings  a and  c along  the  x and  z directions,  respectively. 

In  this  case,  the  height  of  the  lattice  is  taken  to  be  ~ 19^32"- 
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|Ei|/(V0/L) 


Fig.  C-  4 INDUCED  ELECTRIC  FIELD  INTERNAL  TO  DIELECTRIC 
SPHERES  IN  A THREE-DIMENSIONAL  PERIODIC 
LATTICE 
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Fig.  C-5  INDUCED  ELECTRIC  FIELD  INTERNAL  TO  PROLATE 
DIELECTRIC  SPHEROIDS  IN  A THREE-DIMENSIONAL 
PERIODIC  LATTICE 


INDUCED  ELECTRIC  FIELD  INTERNAL  TO  FINITE  CYLINDRICAL 
DIELECTRIC  OBJECTS  IN  A THREE-DIMENSIONAL  PERIODIC 
LATTICE 


* 


Based  on  the  two-dimensional  interaction  constant,  given  by  Equation  (B-ll) 
in  Appendix  B,  the  induced  internal  electric  fields  are  calculated  for  di- 
electric objects  of  prolate  spheroid,  sphere  and  finite  cylinder.  The  re- 
sults are  presented  in  Figures  C-7  to  C-9  for  the  case  when  = 2L/5 

and  %2  = L/5- 

One  lattice  structure  of  relative  importance  is  the  cubic  array  con- 
sisting of  identical  spheres  with  a radius  rQ.  As  a cubic  array,  the  lattice 
element  spacings  are  simply  a = b - c.  Shown  in  Figure  C-10  is  the  induced 
internal  electric  field  for  the  case  when  = £3  = L/4  and  ~ L/2 . The 
induced  electric  field  is  plotted  as  a function  of  the  relative  permittivity 
for  two  values  of  £ (= 

4.  REMARKS 

Simple  expressions  have  been  obtained  for  the  induced  electric  field 
internal  to  the  biological  objects  of  either  a three-dimensional  or  a two- 
dimensional  periodic  lattice.  The  discussion  is  based  on  the  theory  of 
artificial  dielectric.  A limitation  of  this  simple  theory  is  that  the 
object  size  and  the  element  spacing  must  be  small,  being  not  greater  than 
0.1  of  a wavelength.  For  our  present  problem  at  extremely  low  frequencies, 
the  results  should  provide  sufficient  accuracy,  because  the  dimension  of 
the  biological  objects  and  their  element  spacing  are,  indeed,  very  small 
compared  to  a wavelength. 

The  presentation  given  in  this  Appendix  is  based  on  the  analysis  being 
generally  applicable  to  many  problem  areas,  and  not  dependent  upon  the  spe- 
cific details  of  the  simulator,  lattice  structures,  and  the  geometric  shape 
of  the  objects.  Hence,  the  results  are  useful  in  the  study  of  electro- 
magnetic coupling  and  interference,  electromagnetic  biological  effects  and 
biomedical  testings,  slotted  waveguide  arrays  and  artificial  dielectrics. 
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Fig.  C-7  INDUCED  ELECTRIC  FIELD  INTERNAL  TO 

DIELECTRIC  SPHERES  IN  A TWO-DIMENSIONAL 
PERIODIC  LATTICE 
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Fig.  C-8  INDUCED  ELECTRIC  FIELD  INTERNAL  TO  PROLATE 
DIELECTRIC  SPHEROIDS  IN  A TWO-DIMENSIONAL 
PERIODIC  LATTICE 


Left  Scale 


INDUCED  ELECTRIC  FIELD  INTERNAL  TO  FINITE  CYLINDRICAL  DIELECTRIC 
OBJECTS  IN  A TWO-DIMENSIONAL  LATTICE 


|Ej|/(V0/L) 


€2/€0 


Fig.  C- 10  INDUCED  ELECTRIC  FIELD  INTERNAL  TO  DIELECTRIC 
SPHERES  IN  A CUBIC  PERIODIC  LATTICE. 
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